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ABSTRACT 

We discuss a model for the central region of radio-quiet AGN in which the coronal 
accretion is the source of energy for X-ray emission. We consider two-layer geome- 
try and we construct solutions taking into account radiative interaction and pressure 
balance between the disc and the corona. We define three models of the disc-corona 
transition region. Wc use two descriptions of the angular momentum transport: ei- 
ther the angular momentum is transported locally both in the disc and in the corona 
or all the angular momentum is transported within the disc only. We employ the 
a-description of viscosity within the two layers but we discuss also models parame- 
terized by a fraction of energy dissipated in the corona. Both layers are treated using 
vertical averages but with a clear distinction between the values of the pressure and 
other relevant functions at the disc surface and the mean values inside the disc. 

Both descriptions of the angular momentum transport lead to similar conclusions 
about the fraction of energy liberated in the disc although they predict strongly dif- 
ferent properties of the disc interiors. 

The coupled disc/corona system has one family of solutions at low accretion rates, 
two families of solutions at intermediate accretion rates and no solution for high accre- 
tion rate. Of the two families of solutions, the first one has weakly dissipating coronae 
and the second one has coronae which dominate energetically the system. The two 
solutions merge at sub-Eddington accretion rate because a corona in hydrostatic equi- 
librium can dissipate only a limited amount of energy, which depends on the viscosity 
parameter. 

At higher accretion rates an outflow from the corona in the vertical direction 
should be taken into account. 

Key words: accretion, accretion discs - galaxies: active - galaxies: Seyfert - X-rays: 
galaxies 



1 INTRODUCTION 

The overall X-ray spectrum of a typical Seyfert galaxy con- 
sists of the primary featureless component (a power law with 
energy index ~ 0.8 — 1.0) and a component being a reflection 
of the primary emission from cold gas (Pounds et al. 1990 
and subsequent papers; see Nandra & Pounds 1994 for re- 
cent update). Approximately half of the primary radiation 
is seen unprocessed and half of it is reflected. 

Recent data on X-ray and 7-ray spectrum of Seyfert 
galaxy NGC 4151 from OSSE instruments on board of 
Compton GRO suggest that the high frequency emission 
is predominantly of thermal character, since the high en- 
ergy cut-off is observed at only ~ 120 keV (Maisack et al. 
1993) although in the second source measured (IC 4329a) 
the high frequency cut-off of the primary component was 
significantly higher (240 keV < Ec < 900 keV; Madejski et 
al. 1995). This thermal interpretation of the X-ray emis- 
sion Is also strongly supported (although not proved) by the 
shape of the cosmic X-ray background, which is thought pro- 



duced by the overall population of active galactic nuclei (see 
e.g. Zdziarski, Zycki & Krolik 1993). 

Although it is widely believed that accretion onto super- 
massive black hole is the ultimate source of power in AGN 
we do not know how this process proceeds; in particular, we 
do not know where and how the X-ray emission is produced. 
Therefore, a number of scenarios of the generation of X-ray 
emission from accreting black holes have been suggested so 
far: emission from optically thin inner part of accretion disc 
(Shapiro, Lightman & Eardley 1976), emission from shocks 
in the region of jet formation (Henri & Pelletier 1991), emis- 
sion from hot plasma heated by accreting cold blobs (Sivron 
&; Tsuruta 1993) and emission from accretion disc corona 
(Liang & Price 1977). 

Coronal emission is particularly attractive as it is the 
only scenario which clearly predicts that ~ half of the pri- 
mary (coronal) radiation is intercepted by the cool accretion 
disc. Therefore in this paper we concentrate on this general 
scenario. 

Most of the previous results on the accretion disc coro- 
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nae concern three cases: (i) corona heated by an external 
unspecified X-ray source (Begelman, McKee & Shields 1983, 
Ostriker, McKee & Klein 1991) (ii) corona heated by energy 
flux transported from the disc, e.g. through acoustic or mag- 
netohydrodynamical waves (Liang & Price 1977; Bisnovatyi- 
Kogan & Blinnikov 1977; Galeev, Rosner & Vaiana 1979) 
(iii) weak corona heated by viscous forces (Shaviv &; Wehrse 
1986). 

The first approach is interesting for the outer parts of 
an accretion disc but it does not solve the problem of the 
origin of the X-ray emission. The second approach is not 
much promising as only a small fraction of energy can gen- 
erally be transported to the corona (but see, e.g.. Burn &; 
Kuperus 1988) whilst in Seyfert galaxies a large part of the 
energy seems to originate primarily in the form of X-rays 
(Czerny & Zycki 1994; this fraction seems to be smaller in 
bright AGN - quasars - as indicated by an average steeper 
mean Qq/x spectral slope and weaker reflected component, 
e.g. Williams et al. 1992). The third type of models was 
discussed for gas dominated discs (Shaviv & Wehrse 1986, 
Czerny & King 1989) or were based on the assumption that 
the energy generation is proportional to the gas density and 
not to the pressure (Shimura & Takahara 1993). These as- 
sumptions prevented the formation of a strongly dissipative 
corona. 

The fourth option, i.e. generation of all the energy 
within the corona was suggested by Paczyhski (1978). He 
assumed that half of the released energy is emitted in X and 
7 band and that the other half is intercepted by an isother- 
mal disc. The structure of the corona was not discussed in 
the paper; instead, an assumption was made that the disc 
itself is marginally self-gravitating. 

Liang & Thompson (1979) actually calculated the struc- 
ture of the corona. They assumed that only a fraction / of 
the energy is dissipated in the corona; the corona is a two- 
temperature plasma (with Coulomb coupling between elec- 
trons and ions, as in Shapiro, Lightman &; Eardley 1976), 
it is uniform and satisfies the hydrostatic equilibrium. This 
model postulates the equipartition of the radiative and con- 
ductive fluxes, which then uniquely determines the solution. 
Such a corona has a height comparable to the disc itself and 
the density is low enough to keep it optically thin. A more 
general approach to the geometry of the cold and the hot 
gas was discussed by Wandel & Liang (1991). 

Both assumptions (hydrostatic equilibrium and margi- 
nal conductivity) were relaxed by Haardt &; Maraschi (1991; 
hereafter HM91, see also Haardt & Maraschi 1993; hereafter 
1IM93). Instead, they parameterize the solution simply by 
the fraction of energy dissipated in the corona, /, and its 
optical depth, r. Such an approach has the clear advantage 
of being Yery general but it does not allow to predict the 
influence of the existence of the corona on the disc structure, 
and vice-versa. Two recent papers addressed this issue. 

The requirement of hydrostatic equilibrium was reintro- 
duced by Kusunose &; Mineshige (1994). They parameterize 
their solution with the fraction of energy dissipated in the 
corona, /, and calculate the vertically averaged structure of 
the disc and the corona under the assumption that the mean 
pressure is the same in both. However, the mean pressure 
in the corona is significantly lower than the mean pressure 
in the disc (Svenssson & Zdziarski 1994; hereafter SZ94, 
Czerny 1994). This modification allowed for broader range 



of solutions (SZ94) than found by Kusunose & Mineshige 
(1994). 

Nakamura & Osaki (1993; hereafter N093) replaced the 
arbitrary parameter / by assuming that the dissipation in 
the corona can be modeled in the way usually adopted for 
the disc, i.e. that the flux is generated by viscosity which 
scales with the pressure. Also in that paper the hydrostatic 
equilibrium was described assuming the same pressure in the 
disc and in the corona. 

In the present paper we rediscuss the existence and the 
structure of strong coronae within a new frame. We assume 
that coronal accretion itself is the source of energy of the 
corona, whereas the energy released within the main body of 
the disc is subsequently transported to the corona. A similar 
scenario has been recently suggested by Meyer & Meyer- 
Hoffmeister (1994) in the context of cataclysmic variable 
VW Hyi. 

We describe the coronal accretion adopting the popular 
Q-viscosity description. We consider two cases of the angu- 
lar momentum transfer. In the first model (local viscosity 
model - LVM) we assume that the viscosity operates locally 
both in the disc and in the corona and both the disc and 
the corona transport the angular momentum. In the sec- 
ond model (global viscosity model - GVM) we assume there 
exists a global coupling between the two layers so the en- 
tire angular momentum is transported by the disc. We also 
consider models for the spontaneous division of the flow be- 
tween the cold and hot parts. Therefore our models predict 
the ratio of the energy generated in the disc to that gener- 
ated in the corona. This ratio is, in general, a function of 
radius so there is some exchange of angular momentum and 
mass between the disc and the corona. 

Our models are parametrized by the viscosity coeffi- 
cients a in the disc and in the corona, the value of the cen- 
tral mass and the accretion rate. Wo compare our results 
with that of models parameterized by an arbitrarily assumed 
fraction of energy released in the corona, i.e. models which 
arc not l)ase(i on a description of the accretion of corona. 

The plan of the paper is the following. In Section 2 
we give arguments in favour of the coronal accretion model 
and we explain the advantages of using the Q-viscosity de- 
scription of accretion; we discuss the physical aspect of the 
disc/corona transition (Section 2.1) and we formulate a set 
of algebraic equations describing the accretion How in the 
two-layer approximation (Section 2.2-5). In Section 3 we 
discuss the solutions obtained for the LVM. Global viscosity 
model computations are presented in Section 4. Our results 
are discussed in Section -5. 

Preliminary results of this paper were presented at the 
lAU Symposium 159 "Active Galactic Nuclei across the elec- 
tromagnetic spectrum" (Czerny 1994). 

2 LOCAL VISCOSITY MODEL: 
ASSUMPTIONS AND EQUATIONS FOR THE 
TWO VERTICALLY AVERAGED LAYERS 

The physics of accretion is poorly understood as the mi- 
croscopic mechanism of the angular momentum transfer re- 
mains unknown. Most promising, perhaps, is the viscosity 
provided by small scale magnetic field structure which de- 
velops in the disc at the expense of its rotational energy. 
However, computer simulations of this process (Balbus & 
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Hawley 1991, 1992) do not describe usually the saturation 
ol the process and therefore do not give the parameters for 
a stationary situation. 

Therefore, as in the case of stellar convection, we repre- 
sent all the unknown physics by a single parameter a intro- 
duced by Shakura & Sunyaev (1973). This kind of parame- 
terization is very convenient. What is more, some physical 
attempts to describe the global effect of viscosity can be 
reduced to this scaling (e.g. Tout &; Pringle 1992, Canuto, 
Goldman & Hubickyj 1984). 

This assumption was widely used to calculate the ver- 
tically averaged disc structure in AGN (e.g. Ross, Fabian & 
Mineshige 1993, Cannizzo & Reiff 1992, Hure et al 1994a,b; 
see also Frank, King &; Raine 1992). It was also successfully 
used as a local (not vertically averaged) prescription for vis- 
cosity in cataclysmic variables (Meyer & Meyer-Ilollmeister 
1981; Smak 1982; see Canizzo 1994 for recent review) as well 
as AGN (e.g. Lin & Shields 1986; Mineshige & Shields 1990; 
Shimura &; Takahara 1993; Cannizzo 1992). 

Models which are not vertically averaged arc fairly uni- 
form in the vertical direction. However, thermal insta- 
bility of the radiation pressure dominated branch as well 
as the tendency to develop a complex inflow pattern (e.g. 
Rozyczka, Bodenheimer & Bell 1994) may actually lead to 
significant global changes of the disc structure. Therefore, 
it is important to explore the possibility that the resulting 
accretion flow consists of two different layers in the vertical 
direction, and that the heat generation and the radiation 
pressure arc reduced in the optically thick layer, thus possi- 
bly leading to thermally and secularly stable solutions. 

We assume that the inner accreting layer is relatively 
cool and optically thick (a "disc") and that the outer ac- 
creting layer is hot and optically thin (a "corona"). For 
a given accretion rate, a fraction of the accretion flow and 
the corresponding dissipation proceeds through the corona 
and the remaining fraction of accretion flow with the corre- 
sponding dissipation goes through the disc. These fractions 
are predicted by the properties of the two layers and their 
interaction. 

In this section we consider LVM, i.e. we assume that 
the processes mentioned above are parameterized by a vis- 
cosity parameter Qd in the disc and Qc in the corona. In 
this case the disc and the corona act almost independently 
although there is some exchange of angular momentum and 
mass between them, since the fraction of accretion proceed- 
ing through the corona and the disc is a function of the 
distance from the black hole. 

Both layers are in hydrostatic and thermal equilibrium. 
Their interaction occurs through the pressure balance at the 
base of the corona and through the exchange of radiation 
flux. This last interaction is described as in HM91, HM93 
and N093; it strongly influences both layers as the disc pro- 
vides the source of soft photons for the corona whilst its 
structure is influenced by the irradiation, the reduction of 
heat generation and the external pressure. 

We also assume that the local rotational velocity is Ke- 
plerian both in the disc and in the corona. 

The numerical results will be presented only for the ra- 
dius 10-Rsch, which is fairly representative of the innermost 
part of the disc responsible for generating most of the energy 
(Rsch, the SchwarzchUd radius, equals 2 ff.g = 2GAf/c^). 
Approximate analytical solutions at different radii are dis- 



cussed in Appendix B. For the central black hole mass we 
take Af = 10* Mq, and for the viscosity parameters in disc 
and corona = etc = 0.1. We use the dimensionless accre- 
tion rate, 

M 

rn = — : , 

MBdd 

where Mfidd is the critical (Eddington) accretion rate 
JWBdd = — 5 — = , 

assuming the efficiency of accretion rj = 1/12, as it results 
from the Newtonian approximation. 

2.1 Disc/corona transition 

The transition from the cold layer to the hot Compton 
cooled layer is neither arbitrarily positioned nor infinitely 
sharp. As we see from the well studied cases of weak coro- 
nae (e.g. Shaviv & Wehrse 1986, Shimura & Takahara 1993) 
the inversion of the temperature arises because of the drop 
of the atomic cooling efficiency (~ proportional to square 
of the density) at low densities close to the surface whilst 
viscous heating (proportional to the density) is still signif- 
icant. Therefore, close to the surface, the electron temper- 
ature may reach the Compton temperature. The transition 
from atomic cooling to Compton cooling is fairly rapid in 
terms of gas pressure range. 

In principle, full vertical structure of the gas has to be 
solved in order to determine the location of the transition. 
However, some estimates can be obtained in the optically 
thin regime just by studying locally the heating/cooling bal- 
ance of the gas layer under the influence of viscous heating 
and irradiation by the upper and lower gas layers. This 
method is a simple generalization, for the case of (asymp- 
totically) two-temperature medium, of the analysis made by 
Krolik, McKee & Tarter (1981) and Czerny & King (1989). 

As a result of such calculations (sec Appendix A for 
details) we obtain a relation between electron and ion tem- 
peratures vs. the ionization parameter, H, defined here, after 
Krolik et al. (1981), as 



C-fgas 

where F is the radiation flux. The change of gas tempera- 
ture is relatively rapid in terms of pressure range (see Fig. 
1). The solutions along the branch having a negative slope 
are thermally unstable so there should be a switch between 
the lower (cool) and upper (hot) branches. But there is 
no unique way to predict the position of this switch with- 
out computing the structure of the layer and taking self- 
consistently into account the flux varying with the pressure 
(i.e. optical depth). However, identifying this position with 
the upper bend should at least qualitatively well represent 
the properties of the transition layer. 

Therefore we may define the base of the corona (i.e. the 
surface of the disc) by the requirement that 

H = Hmin. model (a) 

It determines uniquely the density at the disc surface if 
the temperature of the surface and the irradiating flux are 
known. We will take this condition as a component of our 
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Figure 1. Relation between the parameter S = -F/cPgas and 
electron (solid curve) and ion (dotted curve) temperatures of the 
plasma. The upper bend on the curve corresponds to S = r-min, 
and determines the conditions at the disc - corona transition in 
our model (a). See Appendix A for details of calculation. 

basic model. 

However one could choose other criteria to define the 
disc/corona interface, and we shall retain two of those, which 
also define uniquely the disc/corona system. 

As our second choice we use the statement done by 
many authors (e.g. N093) that the disc surface density, ps, 
is equal to the disc mean density, pd 



Ps = Pd- 



model (b) 



In the light of the above discussion this condition looks 
rather artificial. 

Another possibility is based on the fact that the opti- 
cally thick part of an o-disc is thermally unstable if dom- 
inated by radiation pressure (T.ightman & Eardley 1974). 
The presence of the corona generally reduces the radiation 
pressure. Therefore we may argue that an initially weak 
corona will be enhanced by the developing disc instability 
until the radiation pressure of the cold layer drops below the 
instability point. Such a disc would be characterized by the 
marginal value of /3 (gas pressure to the total pressure ratio) 
which is a constant resulting from stability analysis, 

l3 = 0.4 model (c) 

(e.g. Shakura &; Sunyaev 1976). 

2.2 Cool disc 

The disc structure is described by the following equations: 
- hydrostatic equilibrium: 



Pup — Pe — ^KPdHd! 

- energy generation: 
3 

Fd = -ad-Pe^K-ffd, 



(2) 
(3) 



Fd 



radiative transfer: 



KpdHd 

equations of state: 



prad ^ , 



PdTe, Te = {iPT^Ia) 



1/4 



prad , 



k 



■psT„ T, = {3P^f/a) 



1/4 



(4) 



(5) 



(6) 



miip, 

and the condition that the radiation flux emitted by the disc 
surface (i.e. at the base of the corona) is equal to the sum 
of the flux generated in the disc and the fraction of the fiux 
generated in the corona and intercepted by the disc: 



3 prad 
4 



Fd + vFcil-a). 



(7) 



Here the symbols have the following meaning: Pup and Pe 
are total pressures at the base of the corona and in the 
equatorial plane, P™^' and P""'^ are radiation pressures at 
the base of the corona and in the equatorial plane, pd, Ps, 
Hd, Fd and Qd are respectively the disc mean density, the 
surface density, the thickness, the flux generated in the disc 
and the viscosity parameter in the disc, whilst Fc is the flux 
generated in the corona and p is mean molecular weight of 
the disc gas. 

The opacity k (the Rosseland mean) as a function of 
density and temperature is taken from Alexander, Johnson 
& Rypma (1983) for logT < 3.8, from Seaton et al. (1994) 
for logT > 4.0, and it is calculated as an interpolation from 
the two tables in the range 3.8 < log T < 4.0. We neglect the 
non-equilibrium influence of coronal X-rays on the opacity. 
This is justified as X-rays arc thermalised close to the disc 
surface at a Thomson depth of the order of unity. The set 
of equations have to be supplementod with the values of r/ 
(the fraction of the coronal fiux directed towards the disc) 
and the X-ray albedo a; we assume i] = 0.5 and a = 0.15 
after HM91. 

Equation (2) reduces to the standard equation of hydro- 
static equilibrium in the disc if the left hand side is negligi- 
ble. It can be neglected if the corona is relatively unimpor- 
tant but for corona-dominated solution the pressure gradient 
within the disc is small as the pressure is mostly provided 
by the weight of the corona and one has Pe « Pup. This 
second limit corresponds to the assumption made by N093 
that the pressure in the disc interior is equal to the pressure 
of the corona; however, in this limit their equation (2) is 
not valid as it is derived under assumption that the surface 
pressure in the disc is negligible. 

Our equation (3) refiecting the local effect of viscos- 
ity is the same as the corresponding equation of N093 al- 
though they claim that the accretion proceeds only in the 
disc. SZ94 adopted different assumptions. A model based 
on an approach similar to their (i.e. Global Viscosity Model) 
is discussed in Section 4. 

2.3 Hot corona 

We assume that the hot corona is a two-temperature 
plasma. A physical description of such a plasma is given 
by Shapiro, Lightman & Eardley (1976). The gas density 
in the corona is taken as constant and the mean pressure in 
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the corona is equal to the pressure at its base. 

Such a corona is described by the following equations 
- hydrostatic equilibrium: 

Fsott 



down - Pc^KHcfn/r 



(8) 



where the last term represents the pressure exerted on the 
corona gas by the soft flux emitted by the disc, 
- energy generation: 

3 



- Compton cooling of electrons: 
Fc = (e« - l)F,oft, 

- equation of state: 



down 



-Pc 



- total corona pressure: 
Ftot 



p, _ pgas 
^down — J^dow 



+ 



(9) 
(10) 

(11) 

(12) 



where the second term on r.h.s. gives the radiation pressure 
on the corona side. It may be written also as [Fsoft + '?(1 + 
a)Fc]/c. 

In equation (10) we replace the amplification factor used 
e.g. by N093, y/{l — y), by a term which is not singular at 
y = 1, but this replacement does not actually influence any 
of the conclusions. 

The Compton parameter y for an optically thin medium 
is given by the formula: 



y = KesPcHc 



4kTe 



(l + ^) 



(13) 



meC^ \ meC^ 

Gravitational energy is converted into kinetic energy of pro- 
tons which cool by Coulomb interactions with electrons, 



3 k 
2 mn 



1 + 



VmeC^ 



1/2' 



where 



2.44 X lO^MnApcTe^^-^ 



InA ~ 20. 



(14) 



(14a) 



The symbols have the following meaning: P^^^^ is the gas 
pressure at the corona base, pc, He, Ti, Te are respectively 
the density, the thickness of the corona, the ion and the 
electron temperatures, jn and fie are the mean molecular 
weight for ions and electrons. 

Equation (8) contains the term fn/r which corrects for pos- 
sible large geometrical thickness of the corona. It is obtained 
by integrating over z the usual hydrostatic equilibrium equa- 
tion: 

dP _ GM z 

dz - r^ + «^V^^TI^^=' 

An approximate form of the corrective term, vaild if Hd <C 



1 + 



-1/2 



(16) 



In our computations we used the exact formula obtained by 
integrating eq. 15 from Hd to He- 



This correction reproduces correctly the behaviour of 
the mean gravitational potential as a function of the corona 
thickness. The potential increases for small values of Hc/r, 
it reaches a maximum at Hc/r ~ 1 and decreases for 
Hc/r 3> 1. This non-monotonic behaviour leads frequently 
to the existence of two solutions for the disc-corona sys- 
tem. However, the solution corresponding to the decreasing 
branch of the potential is dynamically unstable and should 
lead to an outflow. The existence of the maximum gives a 
limit for the amount of energy which can be dissipated. As 
the right hand side of equation (8) saturates for large val- 
ues of Hc/r and the corona is (by definition) optically thin 
(tc = KesPcHc < 1), the energy dissipated in the corona 
cannot be greater than 



r 



n3 2 ; 
-OcilKr /Ke 



which becomes 



m < 0.07 — 



1/2 



(17) 



(18) 



.o.i; Vl07^sch- 

using eq. 19-21 written below. Therefore, for small radii 
and high accretion rates, the corona cannot dissipate an ar- 
bitrarily large fraction of the total energy. It is necessary to 
stress, however, that this behaviour is based on the coronal 
accretion model and that the limit depends on the value of 
the viscosity parameter ac- At the same time it is a strong 
prediction of the model. 

2.4 Relation to global disc parameters 

The structure of the corona depends on the total generated 
flux, Fd + Fc, and on the local Keplerian angular velocity, 
Qk- Both values can be determined at a given radius r 
from the mass of the black hole M and the accretion rate 
M through the usual relations: 



and 



3GMM 



(19) 



(20) 



where f{r) represents the boundary condition at the margi- 
nally stable orbit 

f{r) = 1 - (3/i;sch/r)i/= (21) 

in the Newtonian approximation. 

The disc structure in the two-layer approximation can 
therefore be calculated for given values of the global param- 
eters M, r and M, and of the viscosity coefficients Qd and 
ac. 

2.5 Models parameterized by the fraction of 
energy released in the corona 

In order to better understand the properties which are 
tightly linked to the questionable a-viscosity parameteri- 
zation of the corona, we also calculate models which are not 
based on this assumption. These models are parameterized 
by the ratio ^ = Fd/{Fc + Fd) (^ is simply related to the 
parameter / used by HM91, HM93: ^ = 1 — /; it is more 
convenient if / ~ 1). They are computed without using 
equation (9). 
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3 LOCAL VISCOSITY MODEL: RESULTS 

The most appropriate way to analyze the interaction of the 
corona and the disc is to solve separately the equations for 
the two layers, parameterizing them with ^. Both sets of 
equations give the pressure at the base of the corona. The 
value of ^ is then determined by equaling the pressure pre- 
dicted by the disc equations to that predicted by the corona 
equations, Pup = Pdown- 



3.1 Qualitative discussion 

In model (a) we can obtain an exact expression for the disc 
pressure, Pup (^) by combining eqs. 1, 6 and 7 



(22) 



-a) + ^1 , 

so the values of Pup and of its derivative, dPup/d^, depend 
both only on global parameters. Whether or not there are 
pressure equilibrium solutions depends on the corona pres- 
sure. One may obtain approximate expressions for it in the 
two limiting cases of low ^ and ^ close to 1 (see below). In 
these two limits one can also estimate the disc pressure in 
model (b). Model (c) is more difficult to discuss as it rep- 
resents an intermediate case between the domination of gas 
and radiative pressure in the disc which cannot be approxi- 
mated analytically. 

3.1.1 Weak corona limit 

In the case of a weak corona {(, ~ 1) we can simplify the 
disc equations by assuming that the pressure at the base 
of the corona. Pup, is much lower than the pressure in the 
equatorial plane, Pe, that the disc is dominated by radiation 
pressure and that electron scattering is the only source of 
opacity. 

Under such assumptions the relations between the pres- 
sure, the disc thickness, the density and the surface density 
in the disc dissipating a fraction ^ of the total flux, and in 
the disc without a corona (i.e. ^ = 1; the corresponding 
quantities are marked with SS) are given by: 



Pss, 



Pd=^ pss, 



(23) 
(24) 
(25) 



(26) 



i.e. the disc becomes denser and thinner as the fraction of 
energy released in the disc drops. At a certain point gas 
pressure starts to dominate. Higher density also causes 
the bound-free and bound-bound transitions to become the 
dominant source of the opacity so eq. (23)-(26) are no longer 
valid. 

In model (b) the pressure at the top of the disc, Pup, can 
be derived from eq. (6) by assuming that radiation pressure 
dominates even at the surface: 

Pd+Pc 



so it rises with ^. The values of Pup at ^ = 1 and of its 
derivative, dPup/d^, are therefore both linearly proportional 
to the accretion rate. 

The pressure at the base of the corona calculated from 
the corona equations (8-14), Pdown, goes to zero for ^ ap- 
proaching 1. The equations can be solved analytically in 
this limit but they are complex owing to the difference be- 
tween the electron and ion temperatures. However, combin- 
ing eqs. (8) and (9), introducing the optical depth of the 
corona, Tc = KespcHc and putting fn/r = 1, we can obtain 
an expression for the pressure: 

TcnK(Pd + Pc) V^'^^_^^l/2 (2g) 



Pd, 



3 



which shows clearly the tendency. 

The condition Pup = Pdown (eq. 27 and 28) leads to 
proportionality 



1 - ^ ~ Mac 



(29) 



i.e. with increasing M the corona dissipates an increasing 
fraction of energy until ^ becomes small. 

3.1.2 Strong corona limit 

A detailed analytical computation of the corona dominated 
solutions in model (b) is carried out in Appendix B. Here 
we perform a qualitative discussion in order to understand 
more clearly the behaviour of solutions. If ^ ^ 1 we can sim- 
plify the equations by assuming that the disc is dominated 
by gas pressure and is almost isothermal, i.e. the difference 
between the temperatures at the base of the corona and in 
the equatorial plane is small and of the same order as the 
difference between the pressures (a more precise condition 
for an isothermal disc is ^Td <Si 1, see Appendix B). We reach 
the following expressions showing trends with ^: 

(30) 

(31) 
(32) 



: const. 



Hd = ^PekOk, 



Pe - Pup = ^'(PeKnK)'. (33) 

The exact value of the Pup obtained from this set of equa- 
tions is given by: 

1/2 / , rr. im \ 1/2 

(34) 



P 

up 



(Pd + Pc)fiK 



jQdK 



■ Ps/Pe 



■ PsTs/pdTe 



In model (b) the second term on the right equals 1. In 
model (a) the surface pressure is given by equation (22). 
Therefore in the first case the value of the pressure depends 
significantly on the opacity whilst in the second case it is 
roughly independent of ^ for strong coronae. 

The expression derived from the coronal set of equations 
looks similar in the limit of small ^ and small luminosities 
(in comparison with the limit given by equation 17) 

1/2 

(35) 



Pd, 



TcOK(Pd + Pc) 



{nil - a) + ^[1 - r,{l - a)]}- 



(27) 



Thus in model (b) the two functions Pup(0 and Pdown (0 
vary initially almost parallel to each other and then diverge 
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unless the opacity varies significantly; if the initial value 
Pdovmi^, <S 1) is higher then Pup(^ ^ 1). there will be a 
crossing point, i.e. an equilibrium corona dominated solu- 
tion. 

This is the case for small accretion rates as then k > Kes 
whilst Tc is very low. However, for high accretion rates the 
correction fn/r (see eq. 8 and 16) to the hydrostatic equilib- 
rium becomes increasingly important and Pdown decreases. 
Above a certain value of M the pressure Pdown ^ 1) is 
smaller than Pup(^ ^ 1) and there is no solution for the 
disc-corona system in pressure balance. 



3.2 Numerical solutions for the pressures in the 
disc and in the corona 

Before we show the results of the above relations Pdown (0 
and Pnp(4) we demonstrate in Fig. 2 the importance of two 
effects included in our corona equations (8)-(14): the cor- 
rection for the geometrical thickness of the corona, fn/rj 
and the radiation pressure term, TcPsoft/c (eq. 8). If neither 
of the efl^ects is included, the solution for Pdown (parame- 
terized by 5) exists for all accretions rates and Pdown is a 
monotonic function of m. However, for m ~ 0.1 the corona 
becomes geometrically thick {Hc/r > 1; see also N093), 
even if it dissipates only 10% of the gravitational energy. 
This shows clearly the necessity of taking into account the 
Jlc/r correction. Its inclusion leads to a non- monotonic be- 
haviour of Pdown (fii): as discussed above, the pressure in 
optically thin corona cannot be arbitrarily large. After a 
maximum near m ~ 0.03, depending on ^, the pressure be- 
gins to drop. This results in an increase of He as the corona 
must dissipate a still increasing amount of energy according 
to the Q-prescription (eq. 9). Consequently, the corona also 
becomes optically thick contrary to the initial assumptions 
(and Pdown increases again), but the solutions still exist for- 
mally for all m. However, when we take into account the 
radiation pressure acting on the corona with Tc ~ 1 the so- 
lutions cease to exist above m = 0.1-0.2 (see also Appendix 
B); this means the corona cannot be in hydrostatic equi- 
librium and a possible outflow of the coronal gas must be 
considered (Czerny et al., in preparation) 

At high accretion rates, advection constitutes an im- 
portant way of transporting energy (Narayan & Yi 1994; 
Abramowicz et al. 1995). Rough approximation for the ad- 
vective flux at radius r is: 

F^,^=PsI!l^, (36) 

Assuming that the corona dissipates bulk of the accretion 
power, we can solve the system of equations (8)-(14) inde- 
pendently of the parameters of the disc (cf. Appendix B). 
This leads to the result: 

~ (o/O.l) [—) (^^J . (37) 

This equation shows that the advective flux should be taken 
into account for m larger than a few percent. 

3.S.1 Pressure in the corona 

The relation Pdown vs. (, for the corona is shown in Fig. 3 
for three values of rh. The general behaviour is as expected 
(Pdown ~ (1 — see eq. 28). For high accretion rates 



Figure 2. Corona pressure, fdown, optical depth, tc, and the 
Hc/r ratio as functions of rh for three values of g: 10~* (left col- 
umn), 0.5 (middle column) and 0.99 (right column). In each panel 
the dashed curve is the solution with neither radiation pressure 
in the corona nor the Hc/r correction included (cf. eqs. 8 and 16) 
. The dotted curve shows the solution with the latter correction 
taken into account and the solid curve is the full solution. Note 
the dramatic change of topology of solutions when the corrections 
are included. 

and lower ^ the correction fn/r is increasingly important as 
He becomes larger than r and the pressure drops because 
of the decrease of the gravitational potential for an enor- 
mously expanded corona. The significance of the effect is 
again demonstrated in Fig. 5 which shows the pressure of 
the corona calculated with fn/r = 1- 

The inclusion of the radiation pressure term in the 
corona equation of hydrostatic equilibrium (eq. 8) results in 
the existence of two solutions for Pdown (and the correspond- 
ing values of other corona parameters) for each value of ^ 
(Fig. 2). One of the solutions contradicts our assumption 
that the corona should be optically thin, except for ^ fa 1, 
so we do not show this solution in the figures. 

For higher accretion rates the two solutions approach 
each other. Above m » 0.1 they form a kind of loop in 
the Pdown vs. ^ plane, i.e. they exist over a limited range 
of ^ near unity. The range of ^ shrinks as m increases, 
shifting towards larger values which mean that the solutions 
only exist if the corona dissipates a small amount of energy. 
Finally, for m > 0.2 no solution for a dissipating corona in 
hydrostatic equilibrium is possible. 

3.2.2 Surface pressure of the disc and solution for ^ in 

model (a) 

In model (a) the transition between the disc and the 
corona is determined by the ratio of hard photons pressure 
to gas pressure, i.e. H, such that Comptonization dominates 
as a cooling mechanism. Since radiation pressure at the top 
of the disc depends only on the partition of energy gener- 
ation, i.e. ^ (eq. 7), the total disc pressure at the corona 
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Figure 3. Corona pressure, -Pdown (solid curves), and disc surface pressure. Pup (dotted and dashed curves), for three values of m in the 
three considered models of disc-corona transition (labels a, b and c, see Section 2.1) and both local and global viscosity models (LVM 
and GVM respectively, see Section 4 for definition of GVM). The number in each panel gives the accretion rate. 



base, PupiOi is determined independently of the disc struc- 
ture. The corona equations provide the pressure Pdown (0, 
and the pressure equilibrium condition Pup = Pdown gives 
the solutions valid for any disc structure, provided that the 
constraint P{z = _ffd) = Pup can be satisfied. 

The disc surface pressure, Pup, calculated as a function 
of (, on the basis of criterion (a) is shown in Fig. 3. For low 
accretion rates and strong coronae the predicted pressure of 
the corona is much higher than the pressure in the disc and 
the equilibrium solution exists for ^ very close to 1. How- 
ever, for an increasing accretion rate the height of the corona 
rises and the effect of the factor fn/r becomes increasingly 
important so the pressure in the corona does not grow as 
rapidly as the disc pressure. For accretion rates around 0.04 
(see Fig. 3) both pressures are comparable and there are ac- 
tually two solutions (cf Fig. 6) for the disc-corona system 
in pressure equilibrium. One of the solutions is character- 
ized by large value of Hc/r (cf. Fig. 2) and it belongs to 
the decreasing branch of the gravitational potential. This 
solution is dynamically unstable. Both solutions merge at 
high accretion rates. 

There is no solution for higher accretion rates since the 
disc pressure becomes larger than the corona pressure. In- 
deed, the disc pressure cannot be smaller than the radiation 
pressure determined by m while the corona pressure cannot 
be arbitrarily large due to the drop of gravitational poten- 
tial. 

From the disc equations one can calculate the depen- 
dencies of the disc parameters on ^. For example, in Fig. 4 
we show the mean disc density, pd, and the surface density, 
Ps, as functions of ^, for an accretion rate rh ~ 0.04. The 
surface density varies monotonically and its value is smaller 
than the mean disc density for very strong as well as for very 
weak corona. In the intermediate range ps is larger than pd- 
It will be shown in a forthcoming paper analyzing the full 
vertical structure of the cold layer, that the density inversion 
is a frequent phenomenon in radiation pressure dominated 
discs (Czerny et al., in preparation). In this case, however. 



Figure 4. Surface, ps, and mean, p^, disc densities as functions 
of ^ in model (a), computed for m = 0.042 and Eluiin = 0.05. 
Density inversions occur in some range of ^. Vertical lines indi- 
cate positions of disc-corona solutions in pressure balance for this 
value of m. 

the reason for the inversion may be the assumed value of 
the surface temperature, equal to Ten. Realistic values of Ts 
should be larger because of irradiation by coronal hard fiux 
and, consequently, ps would be lower. 

3.2.3 Surface pressure of the disc and solution for ^ in 
model (h) 

We discuss this model in some detail as it bears most 

similarity with published ones. It can therefore be used to 
show the importance of the modifications to the physical 
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Figure 5. Corona pressure, -Pdown(0 (solid curves), and disc 
pressures, -Pup(5) (dotted curves) for model (b) and the same 
accretion rates as in Fig. 3 (LVM) but showing the influence of 
Hc/r correction to corona hydrostatic equilibrium equation (eq. 
8) and realistic opacities in the disc case (as compared with simple 
Thomson opacity). To illustrate both effects the thick curves 
show solutions without them while the thin curves present full 
solutions (the same as in Fig. 3) for comparison. 

description of the disc layer introduced in this paper. 

The dependence of Pup on ^ calculated from the disc 
equations, with the surface density equal to the mean den- 
sity, is strongly sensitive to the details of the description. 
PupiO shows peculiar "bends" due to the contribution of 
both gas and radiation pressures. Our basic relation PupiO 
for the disc calculated from eq. (2-7), taking into account 
the contribution of bound-free and bound-bound processes 
to the opacities, is shown in Fig. 3. For low ^ the curve 
is flat, as expected (see eq. 30). For ^ close to 1 and high 
accretion rate the curve is rising (see eq. 27). For ^ close to 
1 but low accretion rates the asymptotic behaviour given by 
eq. (27) is not appropriate as the gas pressure is dominating 
at the base of the corona. 

To demonstrate the influence of the bound-free and 
bound-bound processes we also plot the PupiO relation cal- 
culated with the assumption that the opacity in the disc is 
provided only by electron scattering. The result is shown in 
Fig. 5. The shape of the curve is similar to the case with 
full opacities but the low ^ limit of Pup is now higher by 
at least one order of magnitude, depending on m. This re- 
sults from the fact that in tlic disc interior the "real" k is 
much larger than Kes in this limit (cf eq. 34). For ^ near 
1 and m ^ 0.03 the ratio k/kes < 2 therefore Pup is simi- 
lar to that calculated with full opacities. The existence of 
solutions for the coupled disc-corona system - the crossing 
points in Figs. 3 and 5 - is therefore very sensitive to any 
bend in the pressure curve. 

For low accretion rate we obtain one solution for low ^, 
i.e. with most of the energy release proceeding in the corona. 
This solution has the general tendency that ^soi decreases 
with an increase of m; it disappears entirely for accretion 



rates higher than m ~ 0.03. This solution resembles ba- 
sically the N093 gas dominated branch although there are 
significant differences in the description of the disc structure 
in our case and theirs (see Sec. 2.1). The disappearance of 
this brancli for higher accretion rates is closely related to the 
rapid expansion of the corona and the decrease of the mean 
coronal gravitational potential; if this effect is not taken into 
account, i.e. if /a/r = 1 the solution exists even for m = 0.1, 
with ^ always of the order of 10^^ (see Fig. 5, thin dashed 
curve and thick solid curve). 

There are two additional solutions which appear over a 
limited range of m starting from ~ 0.01: one corresponds in 
most cases to ^ 1 (a disc with a very weak corona) and the 
other is corona dominated (^ ~ 0.1). They appear together 
due to the bend in the disc pressure curve and disappear 
also together due to disappearance of the solutions for Pdown 
(Fig. 3), thus creating a loop on ^ vs. rh plane. 

If the bound-free and bound-bound contribution to the 
opacity is neglected, the solutions for the disc and corona in 
pressure balance exist only when rh > 0.03 (Fig. 5) as the 
pressure in the corona is in this case almost always lower 
than the pressure in the disc (see eq. 34 and 35) if we as- 
sume the same value of the viscosity coefficient in the disc 
and in the corona. It would help to increase a in the disc (or 
to decrease a in the corona). However, this would be con- 
trary to what is usually observed in cataclysmic systems, 
i.e. that viscosity in the hot phase is usually higher than the 
viscosity in the cold phase (e.g. Smak 1984); in other types 
of binaries this ratio is close to unity (e.g. Honeycutt, Can- 
nizzo & Robertson 1994). More solutions could also appear 
if we neglect the correction for the height of the corona. 

3.S.4 Surface pressure of the disc and solution for ^ in 
model (c) 

In our model (c) in which we require that the gas to total 
pressure ratio, /3, is independent of the curve Pup(C) is 
very steep (sec Fig. 3). For small values of ^ there is no 
disc solution in hydrostatic equilibrium since Pi^^ and, con- 
sequently Pe, are small due to the /3 = 0.4 condition. On 
the other hand, for ^ « 1 the requirement for /3 implies 
very high values of ps- This is because under the usual as- 
sumption pd = ps the discs are strongly radiation pressure 
dominated. To make the gas pressure comparable to radia- 
tion pressure the mean gas density must be very large since 
the surface radiation pressure is given by eq. (7). 

However, for intermediate ^, similar to the value for 
wliich our model (b) predicts 6 « 0.4, the pressure equilib- 
rium condition is fulfilled and a solution is possible. The 
surface density for this ^ is lower than the mean density. 

3.3 Properties of the disc corona solutions 

It is interesting to study the properties of a disc-like 
accretion flow at a given radius in the log rh vs. log E plane, 
since the slope of the relation indicates the stability or insta- 
bility of the solution (e.g. Meyer & Meyer-Hoffmeister 1981, 
Wandel & fjiang 1991, Kusunose & Zdziarski 1994). Here we 
perform this study at r = 10 -Rsch and we display the values 
of the total surface density (Ed -I- Ec) as a function of rh. In 
tills plot (Fig. 6) we show only the results of computations 
done for models (a) assuming Hmin = 0.05, since they are 
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Figure 6. Accretion rate vs. surface density relation for model 
(a) computed adopting iimin 

— 0.05. The solutions exist only up 
to m 0.085. For comparison, the solid curves show the relation 
in models parameterized by the fraction of energy generated in 
the disc, 5- 

the most realistic (see Section 2.1). For comparison, we also 
show a plot for the disc without corona = 1), otherwise 
calculated exactly in the same way. 

For small accretion rates the model with a corona fol- 
lows closely the branch without corona. The fraction of 
energy generated in the corona is minor and it has almost 
no influence on the disc structure. Above the accretion rate 
~ 0.04 multiple solutions arc possible (see also Fig. 3). They 
form a kind of loop in this diagram. The horizontal branch 
results from the fact that the relations Pdown(0 ^^'^ ^up(^) 



are parallel for ^ ^ 0.1 (Fig. 3) so, when they cross, the 
entire range of ^ is obtained within a very small range of 
m. The other branch of the log m-log E diagram, along 
which the coronal accretion is important, has actually a neg- 
ative slope which usually suggests some kind of instability 
although a more detailed discussion of stability is necessary 
to confirm it (see e.g. N093). 

For higher accretion rates (above m fa 0.085) there is 
no solution for the disc-corona system. The pressure in 
the corona is lower than the pressure required at the disc 
surface for any partition of the dissipation. It is related to 
the fact that the corona being in hydrostatic equilibrium can 
dissipate only a limited amount of energy (see eq. 17). 

Although the solution for the disc-corona system is not 
unique when expressed as a function of accretion rate, the 
fraction of energy liberated in the disc varies monotonically 
along the loop. The same is true for all the other parameters 
except for very low accretion rate: the dependence of the 
temperature of the corona, of its geometrical thickness and 
the Compton parameter on the accretion rate, are shown in 
Fig. 7. The change at low accretion rate is connected with 
the fact that the solutions become single-temperature, i.e. 
the proton temperature is equal to the electron temperature 
is that range. 

All these results depend on the assumed value of the vis- 
cosity parameter a. A decrease of its value lowers the curves 
in the log m-log E diagram, since the maximum accretion 
rate for a disc-corona system with an accreting corona in 
hydrostatic equilibrium depends linearly on a (see eq. 18). 
Lower values of a also lead to a decrease of the optical depth 
and an increase of both the ion and electron temperatures 
in the corona, as we can see from the analytical discussion of 
Appendix B. Eq. (18) and the relations given in Appendix 
B also indicate how the results vary with radius. Although 
for larger radii the limit for the accretion rate becomes less 
stringent the role of advection term increases (see eq. 37). 

Results for model (b) plotted on log rh vs. log E diagram 



Figure 7. Coronal parameters for solutions of model a) shown as functions of the accretion rate, m: 5, electron and ion temperatures 
Te and Ti, comptonization parameter y, and Hc/r ratio. 
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would form a number of branches, including a long corona- 
dominated branch with a positive slope covering a broad 
range of accretion rates. Such a solution resembles the so- 
lution of N093. Although this idea might be attractive the 
existence of this branch is not well justified physically; it 
simply reflects that fact that the surface density in case (b) 
is calculated as a mean density and depends on the amount 
of flux generated in the disc for ^ -C 1, whilst in model (a) 
this surface density is calculated from the ionization param- 
eter and is roughly independent on ^ in this regime. 

Results for model (c) plotted on m-E diagram would 
be positioned along a straight vertical line. This can be 
most easily seen from eqs. 2-5 supplemented with Pg^'^/Pe = 
const condition. Assuming for simplicity Pup ^ Pe and 
^up** ^ Pe^'^ "ii<^ obtains a solution for all disc parameters 
which is independent of Fd and hence of m. 



3.4 Comparison with models ^ = const 

The range of solutions parameterized by ^ is bounded in 
m (Fig. 6). The upper limit reflects the impossibility of 
obtaining solutions in hydrostatic equilibrium even if the a- 
prescription in the corona (eq. 9) is not used. There is also 
a lower limit coming from the requirement of a large optical 
thickness in the disc which is not fulfilled when rh and ^ are 
small. 

Solutions parameterized by constant ^ form a series of 
curves on the m-E diagram whose slope indicates stability. 
One sees (Fig. 6) that weak corona solutions are still unsta- 
ble due to the dominance of radiation pressure (Lightman- 
Eardley instability) in the disc above certain rh. The critical 
value of rh increases with decreasing ^ and, for ^ ^ 0.1, the 
disc is always stable. 

3.5 Role of pair production 

At high coronal temperatures pair production may bo im- 
portant. The estimated value of the compactness parameter 
for our corona-dominated solutions is, however, moderate: 



Ic 



(Tt He 



" r {r/Rsci 

1, C < 1, m < 0.05 at 



(38) 



{Hc/r ~ 1, C < 1, m < 0.05 at r = 10J?sch). Therefore 
we have solved only a posteriori the equation of pair bal- 
ance to check how important they can be. For conditions 
corresponding to the disc-corona solutions (Fig. 7) we have 
computed the optical depth due to pairs and compared it 
with the total optical depth of the corona. We proceeded 
in a similar manner as HM93. The balance equation for 
thermal pairs can be written (Zdziarski 1985, HM93) 



'(1-A), 



(39) 



where Tp and Tc are the "proton" and total optical depths, 
respectively, (it is Tc that results from coronal energy bal- 
ance) and 



^ /wwj I y' /a. 



(40) 



The input spectrum (presumably due to Comptonization) 
is defined as the sum of a power law and a Wien peak and 



Figure 8. The corona total optical depth, Tc, resulting from coro- 
nal energy balance (filled squares), the "proton" optical depth, Tp 
(open circles) and the optical depth due to pairs. Tee, (stars) com- 
puted for the solutions of LVM model a) (cf Fig. 6 & 7) 



A^p and Aw are constants multiplying the two parts. Their 
ratio can be expressed as (Zdziarski 1985) 



Aw 
Ap 



Via) 



r{2a + 3) 



(41) 



where Pr^ is the mean probability of scattering. We calcu- 
late a and Pt^ from formulae given by Zdziarski et al. 1994, 
appropriate for optically thin uniform slab (although we re- 
alize that Hc/r ~ 1 means rather conical geometry). The 
terms /pp, /pw and /ww are pair production rates due to 
photon-photon interactions (we neglect any other pair pro- 
duction mechanisms) while /a is the annihilation rate. Pair 
escape is neglected in our calculations. For all the rates we 
employ expressions given in HM93 (note that the formula 
B3 for /a in that work is misprinted), Svensson (1984) and 
Zdziarski (1985). Finally, the input spectrum is normalized 
to the density of Comptonized photons which can be esti- 
mated (HM93) as 



A, 



ic(l+Tc) 



c{Ec) 27vr{r + Hc)' 



(42) 



where Lc is the energy released in the corona, (-Be) being 
the mean photon energy. 

Fig. 8 displays the values of the three optical depths Tc, 
Tp and the optical depth due to pairs. Tee, for solutions of our 
model a). Under the most favourable conditions Tee does not 
exceed 15% of Tc, i.e. the pairs are relatively unimportant. 
The reason for this is a combination of non-relativistic tem- 
peratures (0 = kTe/meC^ < 0.4) and significant Tc ~ 0.5. 
Under these circumstances the compactness parameter re- 
quired for efficient pair production is > 1000 (see, e.g.. Fig. 
3 in HM93). 
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GLOBAL VISCOSITY MODEL 



4.1 Description and equations 

As we actually do not know how the viscosity mechanism 
operates we should allow for a range of possibilities. If the 
viscosity is connected with generation of magnetic field we 
can cither assume that the corona itself generates its own 
local magnetic field - and therefore creates a local viscosity 
mechanism - or assume that the cool disc is the only source 
of magnetic field. In this latter case the magnetic field of 
the disc would simply penetrate into the corona thus pos- 
sibly extracting angular momentum from the hot gas and 
transporting this angular momentum to the disc. The hot 
corona, loosing its angular momentum, can accrete; the cool 
disc through its own local viscosity mechanism, has to trans- 
port the entire angular momentum whilst participating only 
partially to the accretion flow and to the corresponding dis- 
sipation. 

Such a picture can be translated into equations in a 
very simple way using the description of disc and corona 
structure given in Sections 2.2-3. We only require that the 
entire angular momentum (not a fraction) is transported by 
the cool disc. This requirement, for a Keplerian disc, is 
formally equivalent to replacing _Fd in eq. (3) with Fd + Fc 
as it was done by SZ94, leaving the fraction of the flux to 
be transported by radiation (eq. 4) unchanged. 

We cannot (in principle) describe now the corona using 
the equation (9) since the local viscosity do not operate there 
any more. Instead, we can assume that the radial velocity 
is a constant fraction of the sound velocity. If we call this 
fraction Oc and assume that the corona is in hydrostatic 
equilibrium then equation (9) has to be replaced by: 



= AirracVsPcHc; 



Pdo 



Pc 
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(43) 



For a geometrically thick corona {He ~ r or larger) this new 
equation is in practice equivalent to the previous one, the 
only difference being the factor /(r) representing the inner 
boundary condition for an accretion disc (eq. 21). 

The predictions for the disc structure from the GVM 
differ strongly frorTi the predictions based on the viscosity 
mechanism operating locally. The equations of the cold 
layer structure (Section 2.2) lead to a different asymptotic 
behaviour if the dissipation of energy in the disc is small. 
Instead of a trend qualitatively represented in the LVM by 
eq. 23 - 26 we find in this case: 

1 1/2 

(44) 
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(45) 



so the pressure diverges for small dissipation inside the disc. 
It simply results from the fact that the disc has to trans- 
port all the angular momentum whilst dissipating none. The 

structure of the corona predicted on the basis of the limit for 
the radial velocity is not significantly different from LVM. 

4.2 Results 

The numerical computations confirm the above described 
trend (Fig. 3). In model (b) the pressure expected at the 



Figure 9. Accretion rate vs. surface density relation for global 
viscosity model a) (circles). The curve shows LVM for compari- 
son. 



base of the corona is much higher than for a standard disc. 
The surface pressure in model (a) is defined independently 
of the disc structure so it is the same as in LVM (cf Fig. 
3), while in model (c) Pup(^) differs from that in LVM but 
not significantly. The corona pressure, Pdown, is also consid- 
erably modified (Fig 3, note that the panels for LVM and 
GVM do not show the same accretion rates). For a given 
accretion rate it is lower than in LVM and the solutions exist 
only up to m « 0.07, in comparison with rh « 0.2 in LVM, 
for the adopted values of parameters. 

Disc-corona solutions in model (a) have a character 
very similar to those in LVM: a weak corona is obtained 
for small accretion rates. For m ~ 0.02 (if Hmin = 0.05) the 
asymptotic values Pup{(, 1) and Pdown(C <C 1) are com- 
parable and a corona dominated branch of solution occurs. 
In model (b) the solutions exist only in a narrow range of 
m ^ 0.02 and they are disc-dominated. 

On the m-E diagram (Fig. 9) model (a) solutions form 
two branches of negative slopes. The horizontal branch ex- 
tends towards larger E, instead of smaller E as in LVM. So 
in GVM there is no turn characteristic for a-disc solutions 
occurring when the gas and radiation pressures are compa- 
rable. Model (c) solutions are positioned along a straight 
line of positive slope but the averaged quantities (density, 
pressure) are much larger in the global model than in the 
local one. Density inversions do not occur in this model. 

The opacity is high in GVM. Consequently the temper- 
ature gradient is large even if the disc dissipates only a tiny 
fraction of the energy because of the increase in the total 
optical depth by four orders of magnitude with respect to 
a standard disc which dissipates the same fraction of en- 
ergy but transports only the correspondingly small fraction 
of angular momentum. The mass of such a disc becomes 
considerable if the disc extends to more than 100 Rsch and 
it becomes self-gravitating. 
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5 DISCUSSION AND CONCLUSIONS 

In this paper we made an attempt to describe the disc- 
corona system without specifying arbitrarily the partition of 
the energy generation between the two layers, but instead 
predicting it. We considered two cases. In the first one 
the corona as well as the disc is accreting and transporting 
the angular momentum through a viscosity parameterized 
by a (LVM - local viscosity model). In the second one we 
assumed that the entire angular momentum is transported 
vertically to the disc and radially through the disc (GVM - 
global viscosity model). The partition between the disc and 
the corona in our favourite model (a) (see Section 2) was de- 
termined by the critical value of the ionization parameter H. 
In both cases we assumed that the corona can be described 
as a continuous medium. 

The direct conclusions from this study are the following: 

i) at this stage we cannot favour any of the two models 
of angular momentum transport. Both predict stable 
discs for accretion rates below 0.001 of the Eddington 
limit (when gas pressure dominates) and unstable discs 
between 0.001 and 0.1 (with the presence of a relatively 
weak corona). There is a corona dominated (perhaps) 
stable solution in the LVM model but for a very narrow 
range of accretion rates. 

ii) the present modeling of the corona is not suitable for 
higher accretion rates because of the limit on the flux 
generated in the corona if it is in hydrostatic equilib- 
rium. An outflow from the corona in the vertical di- 
rection is expected. This effect, together with radial 
advection, may change the picture for high accretion 
rates, 

iii) when the parameters of the disc-corona system are cal- 
culated as functions of the fraction of energy liberated 
in the disc, they depends (but not critically) on the 
underlying assumptions, 

iv) for final solutions in pressure balance between the disc 
and the corona the fraction of energy liberated in the 
disc depends strongly on details of the description: the 
adopted law of the disc-corona transition, description 
of opacities within the disc, the correction for corona 
geometrical thickness, 

v) our physically justified solutions do not possess a stable 
corona-dominated solution of the type found by N093, 
covering a broad range of accretion rates because such 
a solution (obtained in case (b) - see Section 2) is based 
on the assumption that the surface density is equal to 
the mean density in the disc. The mean density varies 
with the fraction of energy generated in the disc ^ even if 
it is low (i.e. ^ ^ 1) whilst the surface density calculated 
from the ionization parameter H does not depend on ^ 
in this limit. 

The present study did not give a full picture of accretion 
so it is not possible at this stage to discuss the possible 
advantages of a model of continuous corona over a clumpy 
corona recently proposed by Haardt, Maraschi &; Ghisellini 
(1994). 

Future investigations should therefore go into two par- 
allel directions. First, the theoretical approach has to solve 
the structure of the disc-corona system for high accretion 
rates, and to predict its dependence on the values of the 
parameters. Second, the observational approach should try 



to give an insight into the relation between the amplitudes 
of variations of the UV flux (which in the disc model would 
correspond to disc thermal instabilities, e.g. Siemiginowska 
&; Czerny 1989), the relative level of X-ray and UV emission, 
and the Eddington ratio. 

If the dynamical study of the disc-corona system con- 
firms our present result that there is no solution for large 
Eddington ratios, one should expect strong differences in 
the UV-X spectrum of objects accreting close to or far from 
the Eddington rate. Many studies of the broad line emission 
of AGN have shown that quasars and Seyfert galaxies dif- 
fer not only in their luminosity, but also in their Eddington 
ratio, close to unity for quasars, and at least one order of 
magnitude smaller in Seyfert galaxies. An important obser- 
vational issue would therefore be to clarify the link between 
the luminosity of AGN and their other properties. Prelimi- 
nary results seem to indicate that the so-called "reflection" 
component observed in the X-ray spectrum of Seyfert galax- 
ies and due to reprocessing of X-rays on cold matter, is ab- 
sent in quasars (Williams ct al. 1992). In Seyfert galaxies 
the existence of reprocessing between UV and X-ray pho- 
tons is also ascertained by the absence of time lag between 
them (Clavel et al. 1991). Moreover it is well-known that 
the X-ray to UV flux ratio is much smaller in quasars than in 
Seyfert galaxies. All these facts could indicate that there is 
indeed a fundamental difference between the accretion pro- 
cess in quasars and in Seyfert galaxies, and the presence or 
absence of a strong accreting corona above the disc could be 
a clue for this problem. 
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APPENDIX A 

We can estimate from simple algebraic equations the tem- 
perature of an optically thin layer of gas at a density p close 
to that of the base of the corona. Such a layer is irradiated 
both from the top (by upper corona) and from the bottom 
(cool disc) and it also participates in the energy dissipation 
due to the Keplcrian shear and to the presence of viscosity. 

Here we assume that the viscous heating is proportional 
to the gas pressure of ions, 

Q+,, = |anK— Ti. (Al) 
2 mii 



This energy goes primarily to ions and later on it is trans- 
fered to electrons via Coulomb interactions. The difference 
of temperature between electrons and protons in the sta- 
tionary case is given by: 

Q+sc = |— '^eip(Ti-Te), (A2) 

(for definition of i/ei see eq. 14a). 

Apart from viscous heating, the layer is heated by ab- 
sorption of outgoing cool photons, which we describe by the 
Kramers opacity coefficient. Compton scattering of incom- 
ing and outgoing photons gives a net heating or cooling, 
depending on the sign of the difference between the electron 
temperature and the Compton temperature of the radiation 
bath. The complex atomic processes providing an efficient 
cooling mechanism at 10^ — 10^ K can be approximated by 
a simple analytical formula for the cooling function A(Te) 
(Bufr& McCray 1974). 

Values of Ti and Te resulting from the energy balance 
can be conveniently plotted vs. H (see eq. 1). The plot is not 
very accurate as the approximation of the cooling function 
is not satisfactory for lower temperatures but the position 
of the upper bend on the curve is reasonably correct. The 
characteristic value of E corresponding to this bend depends 
most strongly on the adopted value of the Compton temper- 
ature, ft is equal to 0.5 for Tic = 4 x 10** K and 0.05 for lO'' 
K. As the extension of the observed primary X-ray compo- 
nent of AGN spectra suggest rather high temperature, above 
100 keV, we adopted the second value of H in most of our 
computations. 

APPENDIX B: ANALYTICAL SOLUTION OF 
THE CORONA DOMINATED CASE FOR THE 
LOCAL VISCOSITY MODEL (B) 

To understand qualitatively the behaviour of the solutions 
found with numerical calculations, it is interesting to per- 
form an analytical study, which moreover give their depen- 
dence on various parameters. 

Bl General conditions for the existence of the 
corona 

In the corona dominated case (^ <§; 1) tlie system of equa- 
tions is simplified, since Fc = Ft, where Ft is tlic sum F^+Fc. 
We assume the same value of a in the disc and in the corona. 

To simplify the discussion, we limit it to model b), 
which corresponds to the largest range of m giving a corona 
dominated solution. The disc structure is then described by 
the following equations (to avoid confusion, we sliall keep 
the same numeration as in the main text). 

The hydrostatic equilibrium equation (2) is not modi- 
fied. The energy generation equation writes: 

^Ft = ^aPeQKHi, (B3) 

and the radiative transfer equation: 



c{Pe 



prad\ 



(B4) 



KpdHd 

The equations of state (B5) and (B6) are identical to equa- 
tions (5) and (6), respectively, and equation (7) writes: 

3 



(B7) 
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The corona is described by the following equations: hydro- 
static equilibrium: 

?7(1 - a)Ft 



energy generation: 

Compton cooling of electrons: 
Ft = (e* - - a)Ft 



(B8) 

(B9) 
(BIO) 



The following equations, (11), (12), (13) and (14), are not 
changed. Completed with equations (15), (19). (20) and 
(21), and with the condition Pup = Pdown, tlioy allow to 
determine all parameters as functions of M, M and r. In 
the following we shall use Ms, the mass in 10® Mq, mo.os, 
the dimensionless accretion rate in 0.03MBdd (both values 
being typical of Seyfert nuclei), rio, the radius in lORsch, 
/lo = f (f) / f (rio) and qq.i = a/0.1. Expressed in these 
units Ft and Qk become: 



Ft = 1.4 X lG"mo.o3 /lo M, 



-1 -3 -2 - 

r-io erg cm s 



and 
Ok = 



2.2 X lO'^Mo" 



-3/2 



(B19) 



(B20) 



To simplify we shall also assume that ri{l — a) is constant 
and equal to 0.5 (this is true within 20%), jti=0.5, and 
fie = 1. 

One can first show that there is no corona dom- 
inated solution if Pdown is dominated by radiation 
pressure, i.e. by the term Ft/c. If it is the case, and since 
Tc is of the order unity, equation (B8) shows that Fdown 
dominated by the second term on the right - the dynamic 
radiation pressure - and is therefore negative (unless the 
two terms on the right balance almost exactly, which im- 
plies that Ft, and consequently M, is very small according 
to equation [B9]). This solution has no physical meaning, 
since the corona would be bloated by radiation pressure. 

Assuming then that the first term of the r.li.s. of equa- 
tion (B8), PcOk^c /h/j-j is much larger than the second, the 
equations of the corona can be solved independently of ^, 
and lead to the following solution (we assume Ti 3> Te and 
we neglect the second order term in [B13] and [B14]) 



Pc = 

El 
r 



7.1x10" 



3/4 ,1/4 ,^-1 / . f \- 

^0.1 iH/r ("lo. 03/10) 



1/4 , 



-9/8 



„ „„ -7/12 J.-5/12 / . J. n5/12 -1/8 

0.77 Qq / fu/^ (mo.03/10) ' , 



n 1/6 J--1/6 / . J- \l/6 -1/4 

Tc = 0.56 Qo^i fu/^ (mo.03/10) r-io' , 

rr n n w i n9 -1/6 /rl/6 / . i \-l/6 1/4 

Te = 2.9 X 10 / /jj'/^ (mo.03/10) , 



Ti = 3.2 X 10"^"^ a'^l''^ 



fl/6 
J H/r 



(mo.03/10) r 



6/6 ^-6/4 
10 



(B21) 
(B22) 
(B23) 
(B24) 
(B25) 



These expressions show that the corona satisfies easily three 
of the four required conditions: the ion gas is non relativistic 
(Ti < lO^^K), He is not much larger than r, and Tc is of the 
order unity. None of these conditions depend on the mass of 
the black hole. The fourth condition, that the electron gas 
be also non relativistic {kTe < me(? , i.e. Te < 6 x lO^K) is 



more stringent, and in particular it implies that the radius 
of the corona is limited to about lOORsch with this type of 
solution. A consistent computation should take into account 
the corrective terms in equations (13) and (14). However, 
since the compactness parameter drops with r much faster 
than Te increases (sec cq. 38) the importance of electron- 
positron pairs should not be significant at large distances. 

Equation (B22) shows that the non dimensional ac- 
cretion rate is limited to m less than 0.03, in order 
for the corona to be geometrically thin. This limit is 
the same as that given by the dominance of the advective 
flux on the viscous flux in the corona. 

Finally the pressure of the corona is given by: 



n 1 o in6 1-6/12 j.6/12 i^-l,, 

Fdown = 1.8 X 10 q1 ' /j//^ Mg X 

(mo.03/10) 



7/12 -19/8 
''lO 



(B26) 



One should now check if there is a solution for the disc, 
with Pup = Pdown. It is clear that if Pup = Pdown and if 
gas pressure dominates in the corona. Pup is also dominated 
by gas pressure. This is due to the fact that the isotropic 
radiation pressure at the disc surface and the dynamical 
pressure at the base of the corona are of the same order, 
so when it dominates in the corona, it also dominates in 
the disc. We are therefore led to the conclusion that gas 
pressure dominates Pup, if a solution exists. 

Since Pup should be dominated by gas pressure, one gets 
from equations (B6) and (B7): 



Pd = 2.35 X 10 



-8 -5/12 ,5/12 



0.1 



H/r 



-3/4, 



/ . , -.1/3 -13/8 

(mo.03/10) r-jo 



The condition P|^' > PJp'* writes then: 



(mo. 03/] 



5/12 5/8 ^ -. 



(B27) 



(B28) 



Assuming Qo.i, fn/r and /lo equal to unity, it leads finally 
to: 



m<0.12r-i/', 



(B29) 



which does not give a limit on the accretion rate as stringent 
as the advective flux. 



B2 Solutions for the disc 



Equations (B4) and (B7) give: 
2^ 



F, 
c 



(B30) 



and it dom- 



where Td is the optical thickness of the disc. 

The second term on the r.h.s. is » Pup 
inates if ^Td <€. |. This is the general condition for radia- 
tive heating to dominate viscous heating inside the disc (cf. 
Hure et al. 1994ab). Let us examine the two cases. 



B2.1 ^Td > I or P;"'' > PJ^''. 

We are again facing two possibilities. Either Pg is also dom- 
inated by gas pressure, or it is dominated by radiation pres- 
sure. 
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B2.1.1 Pe is dominated by gas pressure. Then: 



fJ-rriH 



\ ac J 



(B31) 



Since Pe and Pup are both dominated by gas pressure, and 
prad ^ PJp**, one deduces that Pe > Pup, and equation 
(C2) becomes simply: 



Pe = Pdi^d^^K- 



(B2bis) 



From the set of equations (B2bis), (B3) to (B7), (B30) and 
(B31), one gets the solution for the disc: 



-7 in-7 O-3/10 ^2/5 -7/10 , r-7/10^, 

Pd = 7 X 10 B ' ^ ' j' Mg X 

/ • /. \2/6 -33/20 



^ = 1.6 X 10-* S^/* aj'/^* /f/f M-^/«x 



/ ■ J. N-l/6 1/16 

(mo.03/10) r^'o , 



|^=3.4xlO-*Bi/*a;/f Mg-i/'* 



/• ■ f \l/6 1/16 

(mo.03/10 j , 



rd = 10^S=/*a,?/V^/;>8'/'^ 



/ . J. -.1/2 -9/16 

(mo.03/10) ' r-^o : 



(B32) 



(B33) 



(B34) 



(B35) 



where B = k/kss is generally of the order of 10. 

An important condition is that the disc stays optically 
thick, in order to radiate a quasi blackbody continuum in 
the visible-UV range. With this solution, we see that the 
disc is geometrically thin and optically thick. 

The value of ^ is consistent with our basic assumption 
^ ^ 1. It depends very little on the value of the physical 
parameters, and it is of the order of 10^^. Its exact value 
can be determined by a numerical computation giving B. 

Computing the value of Pup from equations (6), (B27) 
and (B32), one gets: 



Pup = 3 X 10' B-^"° a-l'^° M-"^° 



(mo.03/10) r-10 



2/5 -33/20 



(B36) 



which shows that Pup increases with ^. 

Now we have to check the assumptions made to estab- 
lish this solution. 

The condition ^Td ^ I gives: 



2 ^ n 1 ^! d2 -1/24 _f35/24 /• j. xl/3 -1/2 /■do'7\ 

- < 0.16 B Q(, / fjj ' (mo.03/10) ' ' . (B37) 



3 

/ ■ J- \l/3 

g ~ --U.1 JH/r ("iO.03/10) ''10 

For B = 10, and assuming again ao.i = fn/r = /lo = 1, it 
leads to a condition on the accretion rate: 



m>2 10"'' rl(,\ 



(B37bis) 



which is easily satisfied at any radius. 

The condition that Pe is dominated by gas pressure 
writes: 

150 f^y; (mo.o3/io)-'/' rli' > 1, (B38) 

which gives a condition on the accretion rate (again Qo.i, 
fn/r and /lo =1): 



m < 60 rli' 



(B38bis) 



This condition is strongly dependent on the value of B, but 
even for B ~ 10, it gives reasonable accretion rates m <C 
O.SrJg^ for any value of the radius. 

In conclusion the solution exists at all radii, for 
any accretion rate m smaller than 0.03. It corre- 
sponds to ^ ~ IQ-'^ and to Pup increasing with ^. 

B2.1.2 Pe is dominated by radiation pressure. Then 

Pe = (B39) 
c 

Equation (C2bis) is still valid, since again Pe is ^ Pup. From 
the set of equations (C2bis), (C3) to (C8), and (C39), one 
gets the solution for the disc: 

Pd = 1.5 X 10"" C'^ ao.l Mg"^ 



/ ■ r ^-2 3/2 

(mo.03/10) , 



(mo.03/10) r-^o 



7/6^25/16 



(B40) 



(B41) 



^ = 1.6 X 10-^ B-'/' a-l^'' rJ'r Ms'^' 
r ■'^/^ (B42) 

/ . , -,-1/6 9/16 

(mo.03/10) , 



(B43) 



I ■ , -.1/6 -1/16 

(mo.03/10) ■rjo ' . 

The value of ^ is consistent with our basic assumption ^ <C 1 
and the disc is geometrically thin and optically thick. But 
^ is larger than in the previous solution, since it is ~ 10^^ 
and, contrary to the previous case, it depends strongly on 
the radius. So the solution should disappear for rio larger 
than a few units. From equations (6), (B27) and (B40), we 
get: 

Pup = 6.4 X 10' S"^ Qo.l Mg-' 

/ . t x-2 -3/2 

(mo.03/10) rjo' , 

which shows that Pup decreases now with ^. 

2 

3 1 



(B44) 



The condition (ja ^ I, gives: 



2^11 in3 D-1 -1 / • t \-l 3/2 

- < 1.1 10 B Qo.i (mo.03/10) 



(B45) 



For B = 10, ao.l = /h/t- = /lo = 1, it leads to a condition 
on the accretion rate: 



■ ^ c 3/2 
m <C 6 , 



(B45bis) 



which is satisfied at any radius. 

The condition that Pe is dominated by radiation pres- 
sure writes: 



/in ^ D 7/12 j,5/12 / . J. -,7/12 -7 

40 > S ao\ /jj'/^ (mo.03/10) r-io 



7/12 -7/8 



(B46) 



which gives a condition on the accretion rate (with qo.i 

fll/r = /lO = 1): 



m «C 0.4 r,^ B 



3/2 n-12/7 



(B46bis) 



This condition is again strongly dependent on the value of 
B. For B ~ 10, it corresponds to a quite small accretion 



rale, rh <^ 10 ^rfg^. However the limit on the accretion 
rate increases rapidly with the radius. 
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In conclusion the solution exists at all radii for 
small accretion rates (rh ^ 10~^). It should dis- 
appear for rio larger than a few units because ^ 
becomes too large. The solution corresponds to 
^ ~ 10"^ for rio ~ 1 and to Pup decreasing with ^. 



B2.2 ^Td < I or 



2Pt/3c. 



Pip", so the disc is 



This hypothesis implies that Pi'^' 
isothermal, Pe is dominated by gas pressure, and 

Pe = Pup- It is the case considered by N093. It leads 

however to an inconsistency, since one gets i^a = from 
equations (Al) and (A2). This is due to the fact that the 
density at the surface of the disc, ps, is actually not exactly 
equal to the density at the equatorial plane, pd- Let us call 
A = ps/pd < 1- Equation (Al) then writes: 



(B2ter) 



1-A 

Note that N093 avoid the singularity by setting simply Pe = 
^KPdH^. Their equation is an acceptable approximation to 
(B2ter) only if A is small, and it is not consistent with the 
condition Pup > Q^pdH^. 

The set of equations (B2ter), (B3) to (B7) yield the 
solution for the disc: 

-4 ^/^^ ^7/12 ^-1/8 

(B47) 



^ = 2.8 X 10" 



A 



, ■ t \-7/24 1/4 j.6/12 

(mo.03/lo) T-io iair ' 



Pup, Pe being dominated by gas pressure. This solution 
exists at large radii and relatively large accretion rates. The 
completely isothermal solution with Pe = Pup, as the one 
proposed by N093, does not exist. 

(3) At "intermediate" values of ^, ~ 10~^, a solution ex- 
ists for a geometrically thin and optically thick disc with 
Pe > Pup, Pe being dominated by radiation pressure. This 
solution is limited to small accretion rates and to small radii. 

This rather approximate study leads to several conclu- 
sions which agree with those obtained from the numerical 
computation for the same parameters (radius, mass, viscos- 
ity parameter). Obviously it allows only to know if there is 
a solution of the disc-corona model in an asymptotic case, 
when either the radiative, or the gas pressure, is negligi- 
ble. Other solutions might exist in intermediate cases when 
none of the pressures dominates completely, but this can be 
checked only by detailed computations. 



This paper has been produced using the Blackwell Scientific 
Publications TJgX macros. 
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Td = 10^S^/^^ao.?/l=Ml/« 



(mo. 03/10) 



11/24 ^-1/2 ^6/12 
H/r 



' 10 



(B49) 



As Td must satisfy the condition ^ra ^ |, then finally Td 
should verify 1 <C Td The first condition is easily 

satisfied, but the second one gives: 



0.36 B 



-A 



I . , sl/6 -1/4 ^ ^ 

(mo.03/10) ' < -, 



(B50) 



which can be satisfied only if A is very close to unity. We 
conclude that this solution is not consistent, and the 
disc is never completely isothermal. 



B3 Conclusion 

To summarize, we have found that two solutions for the 
corona dominated case, corresponding to two different values 
of ^, exist for small accretion rates, one being limited to 
small radii, and the other to 10^ gravitational radii. For 
larger radii, the corona becomes relativistic and pair effects 
should be taken into account. 

(1) We have shown that the only case where an optically 
and geometrically thin corona can be maintained in hydro- 
static equilibrium is when the radiative pressure at the sur- 
face of the cold disc is negligible, and it corresponds to a 
limit on the accretion rate: m <^ 0.12 r^g^, which does not 
give a limit on the accretion rate as stringent as the advec- 
tive flux. 

(2) At "smaU" values of ^, ~ 10"'', there is a solution 
for a geometrically thin and optically thick disc with Pe > 
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